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ABSTRACT 

We present the results of the simultaneous deep XMM-Newton and Chandra 
observations of the bright Seyfert 1.9 galaxy MCG— 5-23-16, which is thought 
to have one of the best known examples of a relativistically broadened iron Ka 
line. The time averaged spectral analysis shows that the iron K-shell complex is 
best modeled with an unresolved narrow emission component (FWHM < 5000 
km/s, EW ~ 60 eV) plus a broad component. This latter component has 
FWHM ~ 44000 km/s and EW ~ 50 eV. Its profile is well described by an 
emission line originating from an accretion disk viewed with an inclination angle 
~ 40° and with the emission arising from within a few tens of gravitational 
radii of the central black hole. The time-resolved spectral analysis of the 
XMM-Newton EPIC-pn spectrum shows that both the narrow and broad 
components of the Fe K emission line appear to be constant in time within 
the errors. We detected a narrow sporadic absorption line at 7.7 keV which 
appears to be variable on a time-scale of 20 ksec. If associated with Fexxvi 
Lya this absorption is indicative of a possibly variable, high ionization, high 
velocity outflow. The variability of this absorption feature appears to rule out 
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a local {z = 0) origin. The analysis of the XMM-Newton RGS spectrum reveals 
that the soft X-ray emission of MCG-5-23-16 is likely dominated by several 
emission lines superimposed on an unabsorbed scattered power-law continuum. 
The lack of strong Fe L shell emission together with the detection of a strong 
forbidden line in the O vil triplet is consistent with a scenario where the soft X- 
ray emission lines are produced in a plasma photoionized by the nuclear emission. 



Subject headings: galaxies: active - galaxies: individual (MCG-5-23-16) - galax- 
ies: Seyfert - X-rays: galaxies 



Introduction 



One of the key issues in high energy research on Active Galactic Nuclei (AGN) is the 
study of the 6.4 keV iron Ka emission line profile, which can provide fundamental diagnostics 
of the physical and dynamical conditions of AGN central engines. The fluorescent Fe Ka 
emission line is a prominent and ubiquitous feature in the X-ray spectra of AGN and it 
is believed to be produced in the innermost regions of the AGN, such as the Broad Line 
Region (BLR), the circumnuclear obscuring torus and/or the accretion disk. The profile of 
the line itself provides direct information on the region from which it originates. If the Fe 
Ka emission line is produced far from the nucleus, e.g. in the putative torus, its profile is 
expected to be narrow, while if the line originates in the innermost part of the accretion flow 
a broad and asymmetric profile is predicted as a result of the special and general relativistic 



effects such as Doppler sh ifts, gravitational redshift and light bending (see lFabian et al.ll2000 



Reynolds fc Nowakll2003l for a review). In the latter case, the shape of the profile itself could 
be used to derive information on the nature of the black hole and accretion disk system. 

The observations with ASCA {Advanced Satellite for Cos mology and Astroph ysics) of 
relativistically broadened iron Fe Ka emission lines in AGN ( iNandra et al.l 119971 ). and in 
particular the d etection of a broad and skewed profile in the long ASCA observation of 
MCG-06-30-15 (ITanaka et al.lll995l ). were considered the first evidence that at least some 
line emission originates in the inner part of the accretion disk close to the central black hole. 
However, the scenario emerging from XMM-Newton and Chandra observations of AGN ap- 
pears to be more complex. Indeed, these observations have shown that only a handful of 
objects show the presence of the relativistically broadened line, while the narrow Fe emis- 



sion line at 6.4 keV is an ubiquitous feature in ra any type I AGN (see iBianchi et al.l 12004 



Reeves et al.ll2004j : lYaqoob fc PadmanabhanI 120041 ). Furthermore, the broad component ap- 
pears to be in general weaker than what was expected given the initial ASCA results and 
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in some case it may be absent (i.e. NGC 4151, ISchurch et al.ll2003l ). These observations 
have also shown that the interpretation of the Fe profiles can be strongly dependent on the 
modeling of the underlying continuum, which c an be comp licated by the prese nce of complex 
absorption (i.e. NGC 3783.[Reeves et allbooi NGC 3516. lTurner et ahlboosh and reflection 
components (see iReeves et al.l 120071 and references therein). Furthermore, red- and blue- 



shifted Fe absorption lines, associated with the presence of infalling or outflowing matter in 
the proximit y of the black - hole, have been detected in the X-ray spectra of QSOs and Seyfert 
galaxies (see ICappil l2006l and references therein). These absorption and emission features, 
together with the complexity of modeling the underlying continuum, makes the study of the 
Fe line profiles more complex and thus feasible only for the brightest objects. 



In this framework, MCG-5-23-16 represents one of the best and more robust exam- 
ples of a relativistically broadened Fe line. MCG-5-23-16 is a nearby {z = 0.008486) 
Seyfert 1.9 galaxy, with a typical 2-10 keV flux of ~ 8 x 10^^^ erg cm~^ s~^ making 
it one of the X-ray brightest Seyfert galaxies. Previous X-ray observations showed that 
the X-ray emission of MCG-5-23-16 resembles the classical spectrum of a Compton thin 
(i.e. A^H < 10^"^ cm Seyfert 2, with a soft excess below 1 keV and a column density 
iVn ~ 10^2 cm-2 flDewangan et al.ll2003l : iBalestra et al.ll2004f ). Higher energ y observations 
(i.e. above 10 keV) performed with the Rossi X -Ray Timing Explorer { RXTE-.vNeayei et al. 



19981 : iMattson fc Weavei]l20041 ) and BeppoSAX flPerola et al.ll2002l : lRisahtill2002ar ) detected a 
Compton reflection component, that was interpreted as reprocessed emission fr om the distant 



mole c ular t orus. A strong broad Fe Ka line was first detected with ASCA (jWeaver et al. 
19971 . Il998l ) with an equivalent width EW ~ 200 eV, which could be modeled with a broad 
relativistic line profile {i ~ 50°, where i is the disk inclination angle) plus a narrow core 
at 6.4 keV (equivalent width EW ~ 60 eV). The presence of both these components has 
been subsequently confirmed with Chandra and XMM-Newton observations. In particular, 
a narrow iron line component at ~6.4 keV has been clear ly revealed wi th the Chandra ob- 
servation; the intensity of this narrow core (EW ~ 90 eV; WeaverlboOll ) was found to be in 
good agreement with the ASCA results. Meanwhile, two short XMM-Newton observations, 
whose summed exposure time was ^25 ksec, confirnied th e presence of a under lying broad 
component with an EW ~ 100 eV (jPewangan et al.ll2003l : iBalestra et al.ll2004l ). However, 
the relatively short exposure time of the past XMM-Newton observations did not allow the 
above authors to put strong constraints on the geometry of the emission region. Indeed 
neither the inner disk radius nor the inclination angle could be accurately derived using only 
the XMM-Newton data. 

In this paper we present an analysis of the iron K line profile and variability from 
simultaneous deep XMM-Newton (130 ksec) and Chandra (50 ksec) observations of MCG- 
5-23-16; the analysis and results of the simultaneous deep Suzaku (~100 ksec) are described 
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m iReeves et al.l (l2007l ). The XMM-Newton and Chandra observations and data reduction 
are described in section 2. In section 3 we present the modehng of the time-averaged XMM- 
Newton EPIC and RGS spectra and the results of the spectral fits of the simultaneous 
Chandra HETG spectra. In section 4 we report the results obtained with time-resolved 
spectral analysis, aimed to assess the possible variability of the iron K emission line during 
the long observation and to investigate the relation (or lack of) between the Fe emission 
line intensity and the flux of the underlying X-ray continuum. In section 4, we additionally 
discuss the appearance of a sporadic absorption feature at 7.7 keV (rest-frame) that is 
indicative of a possibly variable, high velocity (v ~ 0.1 c) outflow. The results are discussed 
and summarized in section 5. 



2. Observations and data reduction 



In December 2005, MCG-5-23-16 was observed simultaneously with many different X- 
ray observatories: Suzaku, XMM-Newton, Chandra and RXTE; in Table 1 we report the 
log of the different observations. In this paper we concentrate on the XMM-Newton and 
Chandra ob servations, while th e Suzaku and RXTE observations and results are described 
in detail by lReeves et al.l (120071 ). 



2.1. XMM-Newton 

MCG-5-23-16 was observed with XMM-Newton on 2005 December 8 for a total exposure 
time of about 130 ksec (see Table 1). The pn, MOSl and M0S2 cameras had the medium 
filter applied; the MOSl and M0S2 were operating in Small Window mode, while the pn 
was in Large Window mode. The XMM-Newton data have been processed and cleaned using 
the Science Analysis Software (SAS ver 6.5) and analyzed using standard software packages 
(FTOOLS ver. 6.1, XSPEC ver. 11.3). In order to define the threshold to filter for high- 
background time intervals, we extracted the 10-12 keV light curves and filtered out the data 
when the light curve is 2a above its mean. This screening yielded net exposure times (which 
also includes a dead-time correction) of 96 ksec, 101 ksec and 103 ksec for the pn, MOSl, 
and M0S2 respectively. For the scientific analysis we concentrated on the EPIC-pn data, 
which have the highest signal-to-noise, and we used the MOSl and M0S2 data to check for 
consistency. Taking into account the brightness of the source (the 2-10 keV count rates are 
7.9 counts s"^, 2.7 counts s~^ and 2.8 counts s~^ for the pn, MOSl, M0S2, respectively), 
we ran the sas task epatplot to check for possible pile-up and we found that both in the pn 
and MOS detectors the pile-up fraction is below 1%. However, since we have good photon 
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statistics, when analyzing the time averaged pn spectrum, we decided to use only the pattern 
zero data (which correspond to single events), which are better calibrated, and we considered 
the pattern 0-4 (single and doubles) when we extract spectra with lower exposure time for 
the time-variability study. 

The EPIC pn source spectrum was extracted using a circular region of 37" and back- 
ground data were extracted using two circular regions with an identical radius (37") centered 
at ~ 4' from the source. EPIC MOSl and M0S2 data were extracted using a source ex- 
traction region of 27" radius and two background regions with identical size (27") selected 
on the nearby CCDs. Response matrices and ancillary response files at the source position 
have been created using the SAS tasks arfgen and rmfgen. Background subtracted data 
were then binned to have at least 50 counts in each energy bin. 



The Reflection Grating Spectrometer (RGS: lden Herder et al.ll200ll ) data have been reduced 



using the standard SAS task rgsproc and the most recent calibration files; the total exposure 
times are ~ 97 ksec for both RGSl and RGS2. The RGSl and RGS2 spectra were binned 
at the resolution of the instrument (AA ~ O.lA). 



2.2. Chandra 

Chandra observed MCG-5-23-16 with the ACIS-S with two relatively short exposures 
for a total of 50 ks (see Table 1). For this study the Chandra observations were made with 



the High-Energy Transmission Grating (HETG: lMarkert et al.lll994l ) in the focal plane of the 
High Resolution Mirror Assembly. The Chandra HETG consists of two grating assemblies, 
a High- Energy Grating (HEG) and a Medium- Energy Grating (MEG); the HEG affords 
the best spectral resolution in the ~6-7 keV Fe-K band currently available (~39 eV, or 
1860 km s-i FWHM at 6.4 keV). The MEG spectral resolution is only half that of the HEG. 
The HEG also has higher effective area in the Fe K band. The HEG and MEG energy bands 
are ~ 0.9 — 10 keV and ~ 0.4 — 8 keV respectively, though the effective area falls off rapidly 
with energy near both ends of each bandpass. 

The Chandra data were reprocessed with CIAO version 3.2.10 and CALDB version 
3.0.1. Spectral redistribution matrices (rmf files) were made with the CIAO tool mkgrmf 
for each arm (—1 and -|-1) for the first order data of each of the gratings, HEG and MEG. 
Telescope effective area files were made with the CIAO script fullgarf which drives the 
CIAO tool mkgarf . Again, separate files were made for each arm for each grating for the 



http:/ /cxc. harvard.edu/ciao 
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first order. The effective areas were corrected for tlie time-dependent low-energy degradation 
of tlie ACIS CCDs using tlie option available in tlie mkgarf tool in the stated version of 
the CIAO and CALDB distribution. Events were extracted from the —1 and +1 arms of 
the HEG and MEG using strips of width ±3.6 arcseconds in the cross-dispersion direction. 
Lightcurves and spectra were made from these events and the spectral fitting described below 
was performed on first-order spectra combined from the —1 and -|-1 orders (using response 
files combined with appropriate weighting), but keeping the HEG and MEG spectra separate. 
The background was not subtracted as it is negligible in the energy ranges of interest. 
Examination of the image of the entire detector and cross-dispersion profiles confirmed that 
there were no nearby sources contaminating the data. 

In the following, unless otherwise stated, fit parameters are quoted in the rest frame of 
the source and errors are at the 90% c onfidence leve l for o ne interesting parameter (Ax^ = 



2.71). Abundances were set to those of I Wilms et al.l (120001 ). 



3. Spectral analysis 

3.1. The XMM-Newton 0.3-10 keV continuum 

To characterize the X-ray continuum of MCG-5-23-16 we first fitted the 2-10 keV 
pn data with a redshifte d power-law model, modified by Galactic (A^h = 8 x 10^° cm~^; 
Dickev fc LockmanI Eoool ) and local absorption. For this initial fit, we ignored the 5.0-7.5 
keV band, where the Fe Ka emission is expected. This model yielded an acceptable fit 
(xVdof= 1066.3/987) with F ~ 1.65 and A^h ~ 1-3 x 10^^ cm'^ However, the extrapolation 
of this model to the whole 0.5-10 keV band did not provide a good fit (x^/ dof= 6024.2/1794). 
Indeed, it left strong residuals in the soft (E < 1 keV) band and, as expected, at the energy of 
the Fe Ka emission line. Furthermore, it was not clear if the relatively flat photon index was 
intrinsic or indicative of the presence of emission due to Compton reflection. To illustrate 
this, in Fig. 1 we show the ratio of the pn data to an absorbed power law model fltted over 
the 2-10 keV energy band (ignoring the 5.0-7.5 keV band), with the photon index flxed to 
the best flt value F = 1.82, derived with detailed modeling of the XMM-Newton and Suzaku 
broad band spectra (see Section 3.3). 



In order to model the soft X-ray emission, we added to this model a soft power-law 
component absorbed only by the Galactic column density. The photon index of this soft 
component was found to be steep, F = 3.13 ± 0.10, and even at the CCD resolution of the 
pn instrument, the power-law model left line like residuals (see Fig.. 2; black data points). 
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In particular an emission line is required by the data (Ax^ = —37) at 0.92 ±0.02 keV with a 
flux of 1.2 X 10~^ photons cm~^ s~^. This emission line and the steep power-law continuum 
is confirmed by the MOSl and M0S2 da ta (see Fig. 2, red a nd green data points) and also 
by the simultaneous Suzaku observation (IReeves et al.l 120071 ). This model gives a 2-10 keV 
observed flux of ~ 8.2 x 10~^^ erg cm~^ s~^ and an observed luminosity of ~ 1.3 x 10^^ 
erg s~^. 



This model is still too simple to describe the overall emission of MCG-5-23-16, because 
it does not account for the line emission seen in the RGS spectra (see below). However, it 
demonstrates that scattering of the nuclear power law continuum is a plausible explanation 
for the soft X-ray continuum spectrum. 



3.2. The RGS spectra: soft X-ray spectrum dominated by emission lines 

In order to investigate if the apparent steep soft X-ray photon-index could be due to the 
presence of emission lines which are unresolved at the pn CCD resolution, we analyzed the 
RGS data. Indeed, thanks to the long exposure (~ 100 ksec), the RGSl and RGS2 spectra 
have enough statistics to attempt a spectral analysis (a total of ~ 17500 net counts between 
RGSl and RGS2). The first inspection of the RGS data reveals the presence of several soft 
X-ray emission lines as well as the energy cut-off at ~ 1 keV due to the rest frame absorption. 
We then rebinned the RGS spectra in constant wavelength bins at the spectral resolution of 
the instrument (AA ~ 0.1 A) and we used the C-statistic (jCashlll979l ) available in XSPECS 



for the spectral fit, because with this choice of binning, we have some bins with less than 20 
counts. 

We first fitted the RGS spectra with the baseline model obtained with the pn spectrum. 
This model consists of two components: a primary absorbed power law and a scattered soft 
power law component absorbed only by Galactic absorption; both photon indices have been 
fixed to the value found for the primary power law component (F = 1.82, see Section 3.3). 
Overall, this model provides a reasonable description of the RGS continuum, however, line- 
like residuals are present below 1 keV. We then kept this model, hereafter referred as our 
AGN baseline model and we tested two different additional components for the emission 
below 1 keV, which are: (a) a r nulti-temperature thermal emission model with variable 



abundances for different elements (IMewe et al.lll985l ) or (6) several unresolved emission lines. 



The first model represents the emission expected from a coUisionally ionized plasma; the 
latter reproduces the emission due to material photoionized by the AGN. 



^http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/manual/manual. html 
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When modeling with the multi temperature model, we found that the data could be 
fitted with two thermal components with k^Ti = 0.4:4:~^q^I keV and ^3^2 = 0.15lQ;Qg keV 
(AC = 52 with respect to the AGN baseline model, for 9 additional parameters). The flux 
and luminosity of these thermal components are F(o.5_2)keV ~ 9.7 x 10"^'* erg cm~^ s~^ and 
L(o.5-2)keV ~ 2.2 X 10"^° erg s~\ which are consistent with possible X-ray emission from the 
galaxy. The data below 1 keV still required a contribution from the scattered power law 
component of the AGN baseline model, which has F(o.5_2)keV ~ 2.1 x 10"^^ erg cm~^ s~^ 
Allowing the photon index of this soft power law to vary, we found that it was no longer 
unusually steep (F ~ 2.1). Although not well constrained, the abundances for N, O and 
Ne required with this model are found to be low, Z ~ 0.2Zq {Zj^ ~ 0.2 Zq, Zq ~ 0.2 Zq, 
■^Ne ~ 0.4^0). In particular iron is found to be underabundant, with only an upper limit 
of 0.2 Zq. This is due to the lack of Fe L shell emission lines and is at odds with the flux 
measured for Fe Ka emission line (see Section 3.3 and Table 3). Furthermore this value is also 
in contrast with the Fe abundance measured with the neutral iron edge from the reflection 
component in the s imultaneous Suzaku data {Zpe = 0.4 ±0.1 Zq at the 90% confidence level; 



Reeves et all 120071 ). 



We then fitted the RGS spectra by adding to the AGN basehne continuum model several 
unresolved emission lines, fixing both the soft and the hard power-law photon indices to the 
value found for the AGN primary power-law component (F = 1.82, see Section 3.3). To 
account for the excess of counts below 1 keV, five lines are required (AC = 70); their fluxes 
are listed in Table 2, together with their EWs which range from ~10 eV (Nvil) to ~ 46 eV 
(Ovil). The most likely identifications are 2^1 emission lines from H and He-like O, Ne 
and N (see Fig. 3). 

We tentatively allowed the lines' width to vary and found the Ovil Hea, the Ovil ra- 
diative recombination continuum (RRC), and the Neix Hea lines to be marginally resolved, 
while the N vil Lya and the O vill Lya lines are unresolved. Although a quantitative mea- 
sure is beyond the statistics of the present data, from the width of the Ovil RRC feature 
(E> 0.739 keV) we derived an upper limit of fee 7" < 24 eV on the recombining electron 
temperature. It is worth noting that this low value indicates that the s oft X-ray emission 
originates in a pho toionized rather than in a collisionally ionized plasma (iLiedahl &: Paerels 



19961 : lLiedahllll999l ). Taking into account the brightness of the O Vlll and the N vil Lya line. 



this photoionized material should also produce RRC features from these ions. Neither the 
O VIII (E>0.871 keV) or the N vil (E>0.667 keV) RRC are clearly detected in the RGS spec- 
tra. However, the upper limit on the fluxes of both these features are fairly high (5.9 x lO^^ph 
cm~^ s~^ and 3.8 x lO^^ph cm^^ s^^for O VIII and N vil RRCs respectively). In particular at 
the energy of the O Vlll RRC feature the continuum of MCG-5-23-16 shows a steep rise due 
to the emergence of the absorbed power-law component, which also complicates a correct 
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deblending of the Ne ix Hea triplet (see below). 

The O VII Hea and the Ne ix Hea lines are both triplets, which with the present statis- 
tics, can not be resolved into their forbidden and resonance components. However it is worth 
noting that for both the O vil Hea and the Ne ix Hea lines the energy centroids are close to 
the energy of the rest-frame forbidden line (see Table 2). This suggests a major contribution 
from forbidden lines in each of the triplets with respect to the resonance lines. 

In the case of Neix Hea, its line energy is close to the low energy photoelectric cut- 
off present in the MCG-5-23-16 spectrum due to the local absorption (A^h ~ 10^^ cm^^). 
Furthermore the energy is also close to the O Vlll RRC (>0.871 keV) and to the Fexvill- xix 
3d-2p blend (0.853-0.926 keV). We then included in the model three more Gaussian lines 
to account for the O Vlll RRC feature and for the decomposition of the Ne ix Hea triplet^]. 
The fit did not statistically improve, and we cannot derive any quantitative information on 
the ratio between the intensity of the forbidden and resonance line. In the case of the O Vll 
Hea we measured a width of cr = which is probably due to the presence of the 

intercombination and resonance component. In order to confirm that the detected emission 
is dominated by the forbidden line, we then added two lines and kept the line energies 
frozen for the forbidden (0.561 keV), the intercombination (0.569 keV) and the resonance 
(0.574 keV) lines. The measured ratio between the flux of the forbidden and resonance lines 
(for the latter we use the 90% upper limi t) is ^1.6; which again is evidence of a strong 



contribution from a photoionized plasma (jPorquet &: Dubaijl2000l ). Furthermore with this 



model we found that the width of the forbidden line is now unresolved [a < 4.7 eV). Finally 
an inspection of the values of the centroid energies of the detected hues (see Table 2) shows 
that there are no strong shifts between the theoretical and observed values. The measured 
shifts of these lines give a value of ^ 1 eV (which for example for the O vill Lya corresponds 
to a velocity ^ 500 km s^^). 

In order to test whether these lines can be explained with emission from optically thin 
gas photoionized by the AGN, we replaced the un resolved emission lines wit h a grid of 



photoionized emission models generated by xstar ( iBautista fc Kallmad l200ll ). which as- 
sumes a r ~ 2 illuminating continuum and a turbulence velocity of = 100 km s~^. We 
found that the RGS data are well explained with this model with an ionization parameter 
log^ = 1.29!!^Q;j^3. We then let the photon index of the scattered power-law vary and we 
found that, although it is not well constrained, the value is now similar (F = l.drtofs) ^ 
the primary AGN power law emission, in agreement with the scattering hypothesis. Though 



■^The energy of the forbidden (0.905 keV), the intercombination (0.915 keV) and the resonance (0.922 
keV) components were fixed. 
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even with this model the abundances cannot be constrained, iron (Zpc = l-l^os^o) ^'^^ 
underabundant with respect to the value obtained with the detailed fitting of the Fe Ka 



line and the Compton reflection hump detected with the Suzaku observation (IReeves et al. 



20071 ). 



As a final test, we applied this best fit model to the pn, MOSl and M0S2 data. We 
kept the abundances fixed to the values measured with the RGS spectra and we let only 
the normalization and the photon index vary. This model is now a good description of the 
soft spectrum and no strong residuals are present. Finally with this model we found that 
the scattered component has a F = 2.3 ± 0.3 and F(o.5-2)keV ~ 4 x 10~^^ erg cm~^ s"^ . 
This corresponds to 0.5% of the un-absorbed flux of the primary AGN component and the 
emission due to the photoionized gas (L(o.5_2)keV ~ 10^° erg s~^) is 0.1% of the AGN emission. 



To conclude, although from a statistical point of view either the multi-temperature 
thermal emission model or the photoionized plasma model (xstar) give simil ar results, dif- 



feren t diagnostics suggest that, as already seen in o ther Seyfert2 galaxies (IBianchi et al. 



20061 : iGuainazzi fc Bianchil 120071 : llwasawa et al.l 120031 ) the soft X-ray emission is likely due 



to photoionized plasma which could be associated with the Narrow Line Region. The prop- 
erties of the soft X-ray emission of MCG-5-23-16 that favor this interpretation are: the lack 
of a strong Fe-L shell emission (IKallman et al.lll996l ). the detection of a narrow RRC fea- 
ture from O VII and the stronger contribution from the forbidden line in the O VII Hea triplet. 



In order to assess the extent of the soft flux, we used the Chandra observations and 
created an image combining the 0.3-1.0 keV photons from the zero-order of both observa- 
tions. We created a point-source function (PSF) using the MARxH Chandra simulator. We 
then fit a model to the Chandra image consisting of a constant component, to account for 
the background, and two Gaussian components with the centroid positions tied together, to 
account for both nuclear and extranuclear emission. The model was convolved with the PSF 
and then compared to the data using the Cash statistic. Initially the image was fit with 
<^x = <^y, i-e., circular Gaussian models. This gave a = 0.39"l'^o;^5 for the nuclear component 
and a = 1.57"1q3^ extranuclear component. The count rates for the two components 

were 1.4 x 10~^ counts s~^ and 1.2 x 10~^ counts s~^. Note that some fraction of these 
extents is likely due to residual error in the aspect solution, which would effectively increase 
the PSF. Unfortunately there are no other on-axis point sources bright enough in the field 
to check this. Allowing the extra-nuclear component to be elliptical did not improve the 



4http://space.mit.edu/CXC/MARX/ 
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fit significantly, however it reduced C by 7.2 for 2 additional parameters (the additional a 
parameter and the rotation angle), which is significant at the 2a level. It is likely, however, 
that this asymmetry may be due to aspect errors. As a consistency check we fitted the 
zero-order image from the 2000 Chandra observation with the same model. In this case al- 
lowing for ellipticity in either component did not improve the fit significantly. In this fit the 
best-fitting parameters were a = 0.56"1q5'2 for the nuclear component and a = 2.2"~^q'q for 
the extra- nuclear component, which is consistent within the errors. This is indicative that 
around half of the soft X-ray emission is due to the central point like source and half to an 
extended component. Assuming the current cosmology (Hq = 71 km s~^ Mpc^^, Qa=0.73, 
and f2m=0.27) the soft X-ray emission of this latter component originates within ~ 0.7 kpc; 
this value is in agreement with a possible association with the NLR; i ndeed this extensio n 
is comparable to the extension of the [O lll]A5007 derived with HST data (IFerruit et al.ll2000l ). 



3.3. The iron K band 

We then considered the hard X-ray emission of MCG-5-23-16, using the dual power law 
continuum as described above and only the Gaussian emission line at ~ 0.9 keV. In Fig. 4 
we show the residuals left by the absorbed power law model (with F = 1.65) at the energy 
of Fe K band. These residuals clearly reveal the presence of a strong narrow core at the 
expected energy of the Fe Ka (6.4 keV) and broad wings, which extend from ~ 5.7 keV to 
~ 7 keV. The pn data also show a narrow emission line at ~ 7 keV, due to Fe K/3 and a 
drop at 7.1 keV probably due to a refiection edge. The presence of this l atter component 



was already suggested with the previous short XMM-Newton observations (jPewangan et al. 



20031 ). however the short exposure of these observations together with the lack of any simul- 



taneous observation above 10 keV did not allow the authors to pu t stron g constraints on this 



feature. We adopted the best-fit model obtained by lReeves et al.l (120071 ) from the simultane- 
ous Suzaku observation for the underlying continuum in order to derive the Fe Ka emission 
line properties. Indeed Suzaku's broad band energy range (0.4-100 keV) allowed these au- 
thors to measure the amount of Compton refiection and thus better constrain the continuum 
in the 2-10 keV energy band. This model is composed of: a primary absorbed power-law 
component with an exponential cut-off at high energies (> 200 keV ) and a component due 



to re flection from neutral material (the PEXRAV model in XSPEC, iMagdziarz fc Zdziarski 



19951 ). The parameters of this reflection component are: the reflection fraction, which is 
defined by the subtending solid angle of the refiector R = Q/2tt = 1.1; an inclination angle 
i = 45° and abundance Z = OAZq. When fitting this model to the pn data we kept the 
values of Z, R and the cut-off energy fixed, since they cannot be determined using the lower 
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energy band pass of XMM-Newton. After including the reflection component, the residuals 
no longer show the deep edge at 7.1 keV, which is well modeled with the reflection com- 
ponent (see Fig. 5; upper panel). With this model we found that the primary power-law 
component has a photon index F = 1.82 ± 0.01, absorbed by a neutral column density of 
Nn = 1.49 ± 0.01 X 10^2 cm-2. 

To model the Fe line we first added narrow Gaussian lines at the energies of Fe Ka and Fe 
K/3. For this latter line we kept the energy fixed at 7.06 keV and tied its flux to be 12% 
of the Fe Ka flux. This model clearly leaves an excess of counts at the energy of the 6.4 
keV Fe Ka line (see Fig. 5; middle panel), which can be accounted for by including a broad 
Gaussian line or a relativistic disk-line component (see section 3.3.2). 



3.3.1. Chandra Observation of the Narrow Core. 



In order to measure the parameters (i.e. strength and profile) of the broad compo- 
nent we first derived the width and flux of the narrow core using the simultaneous Chandra 
observations. We combined the Chandra ±1 MEG and HEG 1st order spectra of both ob- 
servations. The combined spectra were rebinned at the maximum spectral resolution of the 
instruments (AA = 0.012 A and 0.023 A for HEG and MEG respectively) and the spectral 
fit was minimized with the C-statistic (jCashlll979h . 

We then adopted the Suzaku best fit model for the underlying 2-8 keV continuum. Thanks 
to the high resolution of the MEG and HEG instruments, the Chandra HETG residuals 
clearly reveal the presence of a narrow core at 6.4 keV (see Fig. 6), which is best modeled 
with a Gaussian line at 6.40;o;o? keV and EW~ 80 eV (AC = 43). With this model we 

1400 km 



measured a width of cr = 32_]^g eV, which corresponds to a velocity width of a 

and is thus indicative of a possible origin from the molecular torus. 
Taking into account that the measured width of the Fe Ka narrow core could be due to the 
presence of the broad component, we added a second Gaussian line. Although the fit did 
not statistically improve, we found an EW of ~ 60 eV (flux ~ 5.6 x 10~^ photons cm"^ s~^; 
a ~ 0.4 keV) and ~ 70 eV (flux ~ 5.8 x 10^^ photons cm"^ s~^) for the broad and narrow 
components respectively, which are in good agre ement with the val ues measured with the 
XMM-Newton (see Table 3) and Suzaku spectra (IReeves et al.l 120071 ). With this model the 
narrow core is no longer resolved and we can place only an upper limit on the width of 
cr < 43 eV, which corresponds to a FWHM < 5000 km s^^ (at the 90% confidence level) 
in good agreer aent with the upper limit measured with a previous Chandra observation of 
MCG-5-23-16 JSalestra et al.ll2004[ ). 
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3.3.2. The broad Fe Ime 

We then adopted the Chandra upper hmit on the width of the narrow core for the 
XMM-Newton fits. The broad component was first modeled adding a second Gaussian hne. 
We also add ed a Comp ton shoulder at 6.3 keV, with its normalization set to 20% of the 



Fe Ka flux (IMattl |2002| ) : the fit improved with a Ax^ = 44 for 3 additional parameters 
(X^ = 1932 for 1785 dof). The broad Gaussian component (E = Q.22l:oil keV) has an EW 
of 66^18 eV and a width of cr = 0.42to;}o keV; which corresponds to a velocity of ~ 20000 
km s"^ (FWHM~ 40000 km s^^). 



We then tested a relativistic diskline model (DISKLINE in XSPEC; iFabian et al.lll989l ): 
this code models a line profile from an accretion disk around a Schwarzschild black hole. The 
main parameters of this model are the inner and outer radii of the emitting region on the 
disk, and its inclination. The disk radial emissivity is assumed to be a power law, in the 
form of r~'^. For the fit we fixed the outer radius to be 400Rg (with Rg = GM/c^) and the 
emissivity g = 3. Finally we assumed the line to be from neutral Fe Ka. From a statis- 
tical point of view this fit gives a similar result to the model with a broad Gaussian line 
(x^ = 1928 for 1785 dof), however the high velocity inferred from the width of the Gaussian 
line is indicative that the line must be produced close to the central black hole; i.e. within 
lOORg and thus inside the Broad Line Region. With this model we found that the inner 
radius is i?in = 48^2o-^g5 ^^^1 the inclination angle is i = AV^i^"] while the EWs of the broad 
and narrow components are EWDiskiine = 53l}^ eV and EWn = 64 ± 6 eV respectively. If 
the constraint on the disk emissivity is relaxed and a flatter emissivity is assumed (g = 2), 
then a disk inner radius of 6-Rg is allowed by the present data. The ratio between the data 
and this best fit model is shown in Fig. 5 (lower panel), an absorption line near 8 keV is the 
only residual. Adding an absorption line to our best fit model improved the fit (Ax^ = 20 
for 2 additional parameters, \EW\ ~ 30 eV, E~ 7.9 keV). In Fig. 7 the resulting EPIC- 
pn spectrum and the best model components are shown. The parameters derived for the 
diskline did not change significantly; the main difference is a slight reduction in its EW, 
which became now 46l}3 Taking into account that this absorption feature is indicative 
of the presence of an ionized absorber (see section 4.1 and 5.3), we tested if the presence 
of a more complex absorber could mimic the profile of the detected broad component. We 
found that the inclusion of a two layers of absorption, characterized by a high (log^ = 3.7) 
and low (log^ = 2) ionization level, did not impact the detection of the broad component 
and its parameters. 
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4. Variability of the Iron line and continuum 

During the XMM-Newton observation, the 2-10 keV flux of MCG-5-23-16 varied from 
~ 7 X 10~^^erg cm~^ to ~ 9 x 10~^^ erg cm~^ s~^. In order to investigate the pos- 
sible variability of the line properties and continuum shape, we flrst tested if there was 
any clear difference between the spectra extracted when the source was in relatively higher 
and lower flux states. We extracted spectra using 2-10 keV threshold of < 4.6 counts 

(F(2_iokcV) ~ 7.4 X 10~^^ erg cm"^ s~^ ) and > 5.5 countss"^ (F{2-iokcV) ~ 8.8 x 10"-^^ 
erg cm~^ s^^ ). We then fltted both spectra with the previous best-flt model, replacing the 
diskline with a broad Gaussian. We found no evidence of variability of either the broad 
or the narrow component; indeed their normalizations are consistent within the errors (see 
Table 3). Furthermore, the width of the broad line is constant. In order to confirm that the 
line is not strongly variable, we checked the difference spectrum, obtained by subtracting 
the low from the high state data. The difference spectrum can be modeled with an absorbed 
power law with photon index F = 1.80 ± 0.09; no strong residuals are left in either the soft 
band or in the iron K band. 

As a second test to assess the possible variability of the Fe emission complex, we divided 
the XMM-Newton observation into 5 intervals with a duration of 20 ksec each; for the fit we 
tied all model parameters except the normalization of the primary power-law component. 
This model gives a statistically acceptable fit for all the 5 spectra. Fig. 8 shows the 5.5-8.5 
keV residuals to this model for all five intervals; no strong deviations from the model are 
required at the energy of the iron Ka line. We then allowed the normalizations of the narrow 
and broad components free to vary. Fig. 9 (panel a and panel b) shows the fiuxes of both the 
narrow and broad components. There is no evidence of variability during the present obser- 
vation; furthermore, the fiuxes of both components are consistent within th e errors with the 



values measured in the pr evious XMM-Newton and Chandra observations (jPewangan et al. 



20031 : iBalestra et al.ll2004j ). This lack of variability can be explained by taking into account 
that MCG-5-23-16 is not highly variable on either short and relatively long time-scales; in- 
deed the source has remained at a similar fiux level (7-9 xlO^^^ erg cm~^ s^^ ) for the last 
10 years. 

The most striking feature that appears to be variable during the XMM-Newton observa- 
tion is a possible absorption line at ~ 7.7 keV (E ~ 7.66 keV, observer frame). This feature 
is present in the average EPIC pn spectrum but as shown in Fig. 8, it is strongest in the 
third spectrum (40 ksec after the beginning of the XMM-Newton observation). To illustrate 
this in Fig. 9 (panel c) we compare the intensity of this absorption feature during the five 
intervals. For this purpose we modeled the absorption with an inverted Gaussian line and 
fixed the energy to the best-fit value found with the spectral analysis of the third interval 
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(E ~ 7.7 keV, |/| ~ 3.2 x 10"^ photons cm-^ s-\ EW = 52 ± 15 eV). The hne is clearly 
variable and it appears to be strongest during the interval with the slightly higher 2-10 keV 
flux (see Fig. 8, panel d), while it is barely detected in the other four spectra. 

4.1. A variable absorption feature at 7.7 keV 

4.1.1. Epic-pn background and calibration checks. 

Before attempting any further modeling, we performed several tests to exclude that the 
7.7 keV absorption feature is due to inappropriate background, binning or pattern selection. 
The EPIC-pn background near this energy range presents two instrumental lines due to Cu 
(8.05 keV) and Ni (7.48 keV) Ka emission lines and an inadequate background selection 
could in principle cause spurious absorption features. However several arguments exclude 
this possibility. First of all the net count rate (~ 2.56 counts s~^) of MCG-5-23-16 in the 
5-10 keV band is ~ 300 times greater than the background (~ 8 x 10^'^ counts s^^). Second, 
the sporadic nature of the feature is indicative that the feature can not be an artifact of the 
background or calibration of the EPIC-pn. Finally there was no background flaring activity 
during this time interval. We conclude that the feature is not due to instrumental or external 
background. 

In order to exclude the possibility that the 7.7 keV feature is due to a binning effect 
we rebinned the pn data of the third interval with a constant energy binning of 80 eV. This 
choice corresponds to about half of the energy resolution of the EPIC-pn camera in this 
energy range (FWHM ~ 150 eV at 6.4 keV; see the XMM-Newton Users' Handbook, Ehle 
et al. 200qj). As shown in Fig. 10, the residuals left by the time averaged best fit model 
confirm the presence of the feature, thus excluding the possibility that it is an artifact of the 
choice of the binning. 

To exclude a pattern selection effect we then compared the pn spectra extracted with 
the pattern 0-4 and pattern selection criteria. Though the latter has 30% fewer counts we 
found no significant difference in the absorption line parameters (Ax^ = 26 for 2 dof; E~ 7.7 
keV, EW~ 60 eV). Finally the presence of an absorption feature is confirmed by the MOSl 
and M0S2 spectra extracted in the same time interval. Although, due to the lower photon 
statistic, the significance of the absorption line is lower in these spectra, both the flux and the 
energy of the feature are consistent (|/| = 2.2±1.8 x 10~^ photons cm~^s~^ E =7.4±0.2 keV 



^http: / /xmm. vilspa.esa.es/external/xmm_user_support/documentation/uhb/XMM_UHB. html 
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= 6; see Table 4) with the values found with the pn datslj. Furthermore the presence of 
the absorption feature is confirmed by the simultaneous Suzaku observation, indeed a weak 
absorption feature is present in the time averaged spectrum (see panel c of Fig. 6 in Reeves 
et al. 2007). The absorption line is weaker than in the XMM-Newton observation, which 
could be explained if we take into account the apparent sporadic nature of the feature, with 
a dilution effect due to the longer duration time (~220 ksec) of the Suzaku observation in 
the Earth orbit. Although the absorption line is not well constrained in the time-averaged 
Suzaku spectrum, the energy of the line at 7.8±0.1 keV is coincident with the XMM-Newton 
data, while the flux of the line is weaker |/| = 1.8±0.9 x 10^^ photons cm~^ s^^ (see Table 4). 



4.1.2. Modeling the absorption feature of the third segment. 

We first fitted the absorption feature adding a Gaussian shaped absorption line keeping 
its width fixed to cr = 0.1 keV. The addition of this line improved the fit with a Ax^ = 33 
for 2 additional parameters (x^/(io/=1451/1505). The line energy is 7.72 ± 0.06 keV with 
a |EW| of 52 ± 15 eV (|/| = 3.2 ± 0.9 x 10~^ photons cm-^g-i). Leaving the width of 
the line free did not improve the fit significantly (Ax^ = 5 for one additional parameter). 
With this fit we find a = 0.2 ±0.1 keV, EW=78 ± 29 eV and an energy consistent with the 
previous best fit (E=7.71 ± 0.08 keV). In Fig. 11 we show the confidence contour plot of the 
line parameters (rest frame energy and intensity) with the line width left free to vary. We 
also attempted to fit the absorption feature replacing the Gaussian line with an edge due to 
K-shell absorption from partially ionized iron. This model gives a best-fit energy of 7.331q 22 
keV and an optical depth of r = 0.09 ± 0.03. The fit is statistically acceptable; however, 
it is worse than the Gaussian absorption model (x^/c?o/=1465/1505, which correspond to a 
Ax^ = 14 worse compared to the absorption line) and more importantly, it is unsuccessful 
at modeling the residuals at 7.7 keV. 



^Unfortunately the Chandra observations do not overlap with this segment of the XMM-Newton spectrum. 
However two possible weak absorption features appear to be present at the rest frame energy of about 7.3 keV 
and 7.4 keV (see Fig. 6) suggesting possible variability of the absorber, although the statistical significance 
of these features is low. 
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4-1.3. The significance of the detection of the absorpt 



ion line 



By applying the standard two-parameter F-test to the drop in of 33 for the addition 
of an inverted Gaussian at 7.71 keV, we found a null hy pothesis probabil i ty for adding this 
extra component of ~ 4 x 10^*^. However, as discussed by lProtassov et al. the F-test 

applied in this way could overestimate the true significance of the detected absorption line. 
In particular the F-test does not take into account the number of time bins in which the line 



is sea rched as well as the range of energy where the line might be expected (see lPorquet et al. 



To asses t he significance of the detecti on we then perfor r ned M onte Carlo simulations 
as described in iPorquet et al.l (120041 ) and in iMarkowitz et al.l (120061 ) , for a similar case of a 
blue-shifted Fe Ka absorption line detected in the XMM-Newton observation of IC 4329a. 
We assumed as our null hypothesis model the best fit model with no absorption feature, 
and we simulated 1000 spectra with the photon statistics expected for a 20 ksec exposure. 
Each simulated spectrum was then fitted with the null hypothesis model to obtain a 
value and we systematically searched for an absorption line over the 4-9 keV energy range 
stepping the energy centroid of the Gaussian in increments of 0.1 keV, refitting at each step. 
We then obtained for each simulated spectrum a minimum and created a distribution 
of 1000 simulated values of the Ax^ (compared to the null hypothesis model), which was 
used to construct a cumulative frequency distribution of the Ax^ expected for a blind line 
search in the 4-9 keV range. Not a single fake spectrum had a |Ax^| > 33, thus the inferred 
probability that the null hypothesis model was correct is < 0.1%. Taking into account the 
number of intervals (5) into which the observation had been split, we derive that the line 
detection is significant at > 99.5%. 

Finally, we performed a similar Monte Carlo simulation to test the statistical signifi- 
cance of the absorption feature in the MOS spectra. For simplicity we ran the simulation on 
MO SI only and we found that in this case the significance is only ~ 61%, mainly due to the 
lower S/N of the MOS spectra at this energy. However the fact that the absorption feature 
is detected by the pn, both MOS cameras and Suzaku suggests that the feature is likely real 
and not an artefact. 



4.1.4- The ionized absorber model. 



As already discussed fr om X-ray spectros copic observ ation of several oth er Seyfert 
galaxies (i.e. MCG-6-30-15 jYoung et all boosl : NGC 3783, iReeves et al.l 120041 : Mrk 509, 
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Reeves et al 



2003 



APM 08279+5255, IChartas et all l2002l : PG 1115+080, IChartas et al. 



Dadina et al.ll2005l: Mrk 335lLonginotti et alJl2007l : E1821+643 lYaaoob fc Serlemitsosll2005 
IC 4329a. iMarkowitz et al lbood) and QS Os (i.e. PG 12 11+143, IPounds et alibood: PDS456. 



20031 ) a likely candidate for the 7.7 keV feature is blue-shifted K-shell absorption due to He- 
or H-like iron. In particular if we assume that the line is due to H-like iron (Fe xxvi Lya at 
6.97 keV) the observed blue-shift suggests that the absorber is outflowing with a velocity of 
the order of 0.1c. 

In order to obtain a more physical representation of the absorber, we replaced the Gaus- 
sian line with a model comprised of a grid of phot oionized absorbers generated by the xstar 
photoionization code (iBautista fc Kallmanll200ll ). For the absorber we assumed a one zone 
photoionization model with half solar abundances and a turbulence velocity of 1000 km s~^. 
The free parameters of this model are: the column density (A^h), the outflowing velocity 
of the absorber (fout) and the ionization parameter = L/nr^; where L is the ionizing 
luminosity, n is the electron density, and r the absorber distance). To reproduce the ab- 
sorption feature a column density of ~ 8 x 10^^ cm~^ and an ionization state characterized 
by log,^ = 3.7l'lo;3 erg cm s~^ are required with an outflow velocity of 0.09 ± 0.01 c (~ 30000 
km s^^). A plot of this best-fit model is shown in Fig. 12 (model A). The column density 
is not well constrained and we can place only a lower limit of > 2 x 10^^ cm~^. At this 
ionization level, the Fe K-shell absorption is indeed mainly due to Fe xxvi and is consistent 
with absorption from highly ionized iron outflowing at ~ 0.1 c with respect to systemic. 



A lower velocity outflow could in principle be obtained assuming that the feature is 
due to a Fe xxiv ls-3p line at 7.78 keV. This is illustrated in Figure 11 (model B), which 
illustrates an xstar model with a column density of A^h = lO^^cm"^, an ionization param- 
eter of logE, = 3.0 erg cm s~^ and no velocity shift. The absorption line at 7.8 keV indeed 
corresponds to the above Fexxiv 1 — 3 transition. However, in this ionization regime 
{logC^ ~ 2.5 — 3.0), we would expect to detect absorption trough at 6.5-6.7 keV due to a 
blend of the 1 — > 2 transitions of Fe xvill-xxv; however, this is not observed in the MCG- 
5-23-16 spectrum. In particular, at this lower ionization we would expect to see the strong 
absorption due to the Fe xxiv ls-2p absorption line at 6.67 keV, and for logC, = 3.0 erg 
cm s~^ (as illustrated in the figure) we would also expect to see an even deeper absorption 
feature at 6.7 keV due to the Fe xxv resonant absorption. Furthermore in this ionization 
regime {logC, ~ 2.5 — 3.0), several strong absorption lines from iron L-shell (2 3) transi- 
tions as well as He and H-like Si/S K-shell lines are expected near 1-2 keV, which are not 
observed in either the XMM-Newton Suzaku, or Chandra HETG spectra. Therefore this 
lower velocity solution appears to be ruled out. 
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Finally it is possible to have a low ionization K/3 absorption line without strong Ka 
absorption from the same species. This scenario is shown in Figure 11 (model C) for an ab- 
sorber with a column density of iVn = lO^^cm"^ and an ionization parameter of logC, = 1.5 
erg cm s"^. Indeed, at this ionization state the dominant species is iron less ionized then 
Fe XVII; the L-shell is full and cannot produce the Is — 2p absorption between 6.4-7.0 keV. 
However, there is an absorption line from 7.1-7.2 keV due to 1 — > 3 transitions from Fe less 
ionized then Fe xvil. In this scenario a blueshift of ~ 0.08c would still be required to model 
the absorption line at 7.7 keV in the spectrum. Furthermore the column density required to 
model the EW of the K/3 absorption feature is A^h >> 10^^ cm~^, which would introduce too 
much continuum (bound-free) absorption below 6 keV, inconsistent with the pn data. 
Therefore we conclude that the fast (0.1c) high ionization outflow is the most likely model 
to account for the absorption feature at 7.7keV. Moreover, when we compare the 2-10 kcV 
continuum and the neutral A^h measured during the third interval with the other inter- 
vals and with the average spectrum, wc do not find any statistically significant difference 
{ANh < 10^^ cm^^), which rules out the possible presence of variable neutral and/or low 
ionization absorber. 



5. Discussion and Conclusions 

We have presented the results from XMM-Newton and Chandra observations of MCG-5- 
23-16, which are part of a simultaneous campaign conducted in December 2005 also compris- 
ing Suzaku and RXTE observations. The 0.5-10 keV continuum of MCG-5-23-16 resembles 
at first order the canonical X-ray emission expected from a Compton thin Seyfert 2 galaxy: 
an absorbed (A^h = 1.5 x 10^^ cm~^) power law component (F = 1.82), which emerges at 
energy ^ 1 keV, and a steep soft excess, which is well fitted by a power law component plus 
several emission lines from O, Ne, and N. The XMM-Newton observation of MCG-5-23-16 
confirms the presence of the Fe Ka emission complex, which is well described by a narrow Fe 
Ka emission line superimposed on a relativistically broadened component. The simultaneous 
Suzaku observation provided us with an accurate description of the underlying continuum 
which allowed us to perform detailed modeling of the Fe K emission fine complex. Finally 
due to the sufficiently long duration, the XMM-Newton RGS spectra have enough photon 
statistics to investigate the origin of the soft X-ray emission. 
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5.1. The origin of the soft X-ray emission 



The analysis of the EPIC-pn and MOS spectra of MCG-5-23-16 revealed the presence, 
below 1 keV, of a soft excess with respect to the primary nuclear emission. This soft excess 
can be well fitted by adding an un-absorbed power law component to the primary AGN 
emission, the photon index of this power law is found to be steeper (F ~ 3.1) than the 
primary AGN component (F ~ 1-8), and even at the EPIC CCD resolution an emission line 
is detected around 0.9 keV. 



A soft excess below 1 keV is not unusual in obscured Seyfert galaxies like MCG-5-23-16 



see 



Bianchi et al.l l2006l and references therein) and it has been already suggested that it 



could be due to a superimposition of scattered emission into the line of sight by ionized 



gas plus several emission lines : 


rom highly ionized 


He and H like 


elements (i.e Mrk 3 


Sako et al. 


200d; 


Pounds &; Paee 


2005; 


Bianchi et al. 


2005 


: Circinus: 


Sambruna et al. 


2001 


NGC 1068: 


Kinkhabwala et al. 


2OO2I 


Oele et al. 


2003: 


Brinkman et al. 


2002; NGC 4507 



Matt et al.ll2004l ). Key diagnostics to understand the origin of this X-ray emission when a 
high resolution spectrum is available are: the detection of RRC transitions, the detection 
of enhanced K-shell emission lines from H-like and He-like ions, the detection of Fe L-shell 
emission and the ratio between the forbidden and the recombination transition in the He-like 
triplets. 

This long XMM-Newton observation has provided for the first time RGS spectra for 
MCG-5-23-16 with sufficient photon statistics to perform a detailed modeling of the soft 
X-ray emission, allowing for the first time in this object the detection of the O Vll RRC. 
The width of this emission line indicates that the recombining electron temperature is a few 
eV {k-QT < 24eV). T his is suggestive that th e emitting plasma is photoionized rather then 
collisionally ionized (ILiedahl fc Paereld Il996l ) . This is also indicative that the soft X-ray 
emission is probably dominated by scattering of the primary AGN emission rather than due 
to emission from hot gas in the host galaxy , e.g. from starburst activity. Our analysis of 
the RGS spectra of MCG-5-23-16 confi r ms th e detection of the O vil and O vill lines previ- 
ously reported by iGuainazzi fc Bianchi! (120071 ) and the line fluxes are in agreement with the 
measurement obtained with the analysis of these previous short XMM-Newton observations 
of MCG-5-23-16. We cannot exclude on a statistical ground the possible presence of emis- 
sion due to a collisionally ionized plasma. Indeed the spectra can be equally modeled by 
replacing these emission lines with a multi temperature thermal model; which represents the 
emission due to a collisionally ionized plasma. However the "AGN" model (scattered power 
law component plus several photoionized emission lines) is preferred because of the low ion 
abundance obtained {Z < 0.2Zq) in the thermal model, which reflects the lack of a strong 
Fe L-shell emission with respect to the fluxes of the Oxygen lines. 
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We therefore conclude that the most hkely the origin of the soft X-ray emission is due to a 
plasma photoionized by the AGN. This plasma must be located outside the Compton-thin 
absorber and, as already suggested for other Seyfert 2 galaxies, it could be coincident with 
the NLR. 



5.2. The Fe Kct emission complex 



This deep XMM-Newton observation of MCG-5-23-16 confirms the presence of broad 
and narrow iron Ka emission lines, which were reported since the first ASCA observation 
(j Weaver et al.l 119971 ) . The ChandraH.ETG spectrum clearly reveals a narrow line at Ej^j = 
6.40 ± 0.02 keV with a FWHM < 5000 km s"^ and a flux of 5.6 ± 0.7 x 10"^ photons cm'^ 
s^^. The intensity of this component is found to be constant, within the errors, during 
this observation a nd also when compa ring with previous observation (< /^v >= 4.5 x 10 
photo ns cm" 



-2 S-\ 



Balestra et al.ll2004l : /at = 6. 5 ±2.7 x 10"^ photons cm 



-2 s-i. 



Weaver et al. 



19971 ). The constancy in flux of this line together with the limits on the width obtained with 
Chandra are suggestive of an origin from distant matter such as the putative torus. Indeed 
the upper limit on the FWHM corresponds to a distance from the central black hole greater 
then WRg. 

The presence of two Compton-thick X-ray reprocessors respon sible for the two co m- 
ponents of the iron line, suggested since the ASCA observation (I Weaver et al.l 119971 ) is 
confirmed with this deep XMM-Newton observation and with the deep Suzaku observation 
(IReeves et al.l 120071 ). T he geometry inferred for MCG-5-23-16 with this latter observation 
is discussed in detail in iReeves et al.l (120071 ). To summarize, one plausible scenario is that 
we are seeing MCG-5-23-16 through the Compton thin edge of the putative torus, which is 
Compton thick at the plane of the accretion disk. This is in agreement with the inclination 
{i = 4ll^2°)0 of the accretion disk derived by modeling the broad line component with a rela- 
tivistic line profile. The column density of 1.5 x 10^^ cm~^ measured using low energy cut-off 
is thus associated with the thinner absorbing material, e.g. encountered viewing through 
the edge of the torus. We found no evidence of variability of the column density of this ab- 
sorber within this long observation, and also no strong variability is found when comparing 
the column densities measured with the previous observations performed with BeppoSAX 



'^Modeling the combined Suzaku and XMM-Newton spectra, IReeves et al.l (|2007() found an inclination 



of i = 53tl ° which is in agreement within the errors with the value obtained from the analysis of the 
XMM-Newton spectrum. 
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20031: 



flRisalitill2002ar). ASCA f|Weaver et alJll997h . XMM-Newton and Chandra flDewangan et al. 



Balestra et al.ll2004l ). This result implies that the absorber is probably far from the 



central black hol e and there is no ev idence that this absorber is clumpy as suggested for 
other Seyfert 2s jRisaliti et al.ll2002b[ ). 



This deep XMM-Newton observation confirms the presence of a relativistically broad- 
ened iron Ka line; the width derived from modeling this component with a Gaussian profile 
corresponds to a FWHM ~ 40000 km s~^ and is suggestive of an origin in the accretion disk. 
The profile of this component is nearly symmetric and can be modeled, equally well with a 
broad Gaussian or a relativistic profile; in the latter case the derived inner radius is about 
20-4:0Rg. Since the advent of XMM-Newton and Chandra one of the most debated issues 
in the study of the broad iron Ka l ines has been the fraction of AGNs which clearly show 
the presence of a broad component (INandra et al.ll2006l : iGuainazzi et al.ll2006l ). Several au- 
thors have discussed the robust ness of the d etection in some ob jects of a broad iron line (i.e. 
NGC 3516. iTurner et al.lboosl : NGC 3783. iReeves et al.ll2004} ). This controversy emerged 
when observations characterized by high photon statist ics showed the ambiguity of modeling 
the iron Ka line when complex absorption is present (jPounds et al.ll2003l . l2004l ). Indeed, a 
high column density warm absorber can produce curvature in the spectrum at the energy of 
the iron Ka line that mimics the profile of a relativistically broadened emission line. The 
detection of the absorption feature at ~7.8 keV shows that a high column density variable 
absorber (a high velocity, highly ionized outflow) is also present in MCG-5-23-16; this could 
in principle give rise to ambiguity in the interpretation of the broad component. However, in 
the case of MCG-5-23-16 the availability of a simultaneous observation with Suzaku above 
10 keV allowed us to tightly constrain the underlying X-ray continuum and to rule out the in- 
terpretation of the broad line as due to unmodeled complex absorption. Indeed the residuals 
left at the energy of the iron Ka, when we take into account the amount of reflection detected 
with Suzaku, cannot be explained by the effect of complex absorption. Furthermore, once 
the absorber responsible for the feature detected at ~ 7.8 keV is accounted for, either in the 
time-averaged spectrum or in the third segment of the XMM-Newton observation, a broad 
line is still required by the data with a similar EW and FWHM. Note that the ionization 
parameter of the absorber is required to be high and does not introduce additional spectral 
curvature below 6 keV, which hence does not impact the broad iron line modeling. 

The remaining open questions on the origin of the broad line in MCG-5-23-16 are the 
relatively large inner radius derived for the accretion disk, and its lack of variability. The 
former can be explained with several scenarios: the disk could be truncated or missing below 
20Rg, or the inner part of the accretion disk could be so highly ioni z ed tha t the iron is fully 
ionized. However it worth noting that, as shown by iReeves et al.l (120071 ). assuming a flat 
emissivity (g = 2) a inner radius of about 6Rg (in the case of a Schwarzschild black hole) 
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cannot be statistically ruled out. 



The second open issue is the lack of variability of the iron emission line both on short 
and long time-scales. Indeed the flux of the broad component is found to be consistent with 
being constant when using the short-term time-resolved spectroscopy performed within this 
deep observation, and when comparing our result with the long-term flux history presented in 
Balestra et al.l (120041 ). The strength of the broad component appears to be lower during this 
obs ervation with respec t to the value reported since the first ASCA observation (EW~ 200 
eV; IWeaver et al.l 119971 ). however when we take into account the larger errors on the early 
ASCA measurements we cannot exclude the line being constant. 

In MCG -5-23-16 this lack of variability of the iron emis sion line is not so striking as in MCG- 
6-30-15 flMiniutti et al.l 120071 : IVaughan &: Fabianll2004l ). due to the low level of variability of 
the intrinsic continuum (30-40% compared to a factor of 2-3 in case of MCG-6-30-15). 



5.3. The blue-shifted absorption line: a possible high velocity, ionized variable 

outflow 



Perhaps the most interesting result of this long XMM-Newton observation has been the 
discovery of a possibly variable absorption line from ionized iron. The feature appears to 
be transient with a time scale of about 20 ksec and it is detected at an observed energy 
of about 7.66 keV (corresponding to a rest-frame energy of 7.72 keV). As shown the most 
plausible association of this feature is with K-shell absorption from H-li ke iron, which is blue- 
shifte d by ~ 0.1c. Indeed modeling this absorption feature with xstar ( iBautista &: Kallman 



200ll ) requires a column density of about 8 x 10^^ cm~^ and a high ionization parameter 

which implies that the absorption is due to a blueshifted 1 
6.97 keV). The velocity of the absorbing material is found to 



(log^ = 3.7 ± 0.3 erg cm s ^] 
2 transition of Fe xxvi {E = 
hev = (0.09±0.01)c. 



In the last few years, red- and blue-shifted absorption lines associated with the pres- 
ence of highly ionized gas in- and/o r out-flowing at relativistic velocities have been re- 



porte d for severa l AGN (E182 1-H643 lYaaoob fc Serlemitsod l2005l : Mkn 509 IPadina et al. 



20051 : NGC 3516 iNandra et al.l Il999l : iTurner et al.ll2005l : Mrk 335, iLonrinotti et al. 



These ab sorption lines are foun d both in Se yfert galaxies (NGC 3783, iReeves et al 



IC4392a, Markowitz et al. 2006; NGC 1365 



and BAL QSOs (APM 08279+5255. IChartas et al 



Risaliti et al. 



2007h as well as in quasars (PC 121 1+143. IPounds et al.l l2003l: PDS 45 6.lReeves et al.ll2003h 



2002 



2005; Ark 564, 



PC 11 15+080. IChartas et al. 



2007). 



2004 



Papadakis et al. 



20031). 



These absorption systems can also be variable on different time scales, in their ionization 
state and column density, with the most extreme cases being NGC 1365 (iRisaliti et al.ll2005l ) 
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and Mkn 509 (iDadina et al.ll2005l ). From the analysis of the different intervals in which we 
split the observation we can infer that it is unlikely that the variability of the absorber in 
MCG-5-23-16 is due to a change in the ionization state of the outflowing material, otherwise 
we would detect signatures of this absorber in all the spectral slices. A more likely scenario 
is a change in column density of this absorber. 

The present data suggest we are seeing a transient absorber, which could be associ- 
ated with a cloud which sporadically obscures the central source. This "cloud" could be 
the signature of a clumpy absorber located close to the central X-ray source or of matter 
ejected sporadically. Different models have been proposed to explain the powerful outflows 
detected with the recent XMM-Newton and Chandra observations; in particular transient 
red an d blue shifted absorption lines are predicte d in several theoretic al models of failed dis k 
winds (jProga et al.l 120001 : iKing fc Pounds! l2003l ) or an aborted jet f lGhisellini et al.l 120041 ). 
The picture emerging is that these features can provide a direct probe of the dynamics and 
kinematics of the innermost central regions of AGNs. For the absorber detected in MCG-5- 
23-16 the data suggest that this feature appears when the source reaches a relative maximum 
in the intrinsic 2-10 keV flux. However, monitoring the spectrum on longer timescales would 
be required to determine whether there is a statistically firm correlation between the pres- 
ence of the absorber and the source brightness and determine if any duty cicle is present. 
Thus both a scenario where a clumpy absorber or a variable or failed outflow or jet are at 
present possible. Finally, it is worth noticing that the rapid variability of this absorption 
feature is indicative of a compact absorber and rules out a possible z = origin, e.g. due to 
the warm intergalact ic medium (WHIM) o r a loc al hot bubble as claimed along the line of 
sight to other AGN jMcKernan et al.ll2004j . mOM)- 



Before deriving an estimate of the location, mass and energetics associated with the ab- 
sorber we performed a consistency check for the A^h value measured with the xstar model 
and the EW (~ 50 eV) of the absorption line measured with the Gaussiari component. Fol- 
lowing the curve of growth for H-like iron (see Fig. 4 of iRisaliti et al.ll2005l ) we inferred that 
the detected EW requires a turbulence velocity greater than 500 km a~^; a lower turbulence 
velocity would imply a Compton thick absorber and a similar EW of the H-like K/3 line (due 
to saturation of the Ka line), which is not detected (|-EW^| <15 eV at the 90% confidence 
level). On the other hand a turbulence velocity greater than 3000 km s^^ would produce 
a broad absorption feature which would be resolved even at the EPIC-pn resolution at 8 
keV (~ 170 eV). We therefore conclude that the observed EW and line width are in broad 
agreement with the XSTAR estimate of a column density of about 10^^ cm~^, for a turbulence 
velocity cr = 1000 km s~^. Using this value for the column density we can now estimate the 
maximum distance of this cloud or blob from the central black hole using the relation be- 
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tween the ionization parameter, the density of the absorber and the illuminating continuum 
luminosity: Lj^ = nR^; where L is the intrinsic 2-100 keV X-ray luminosity (5.4 x 10^^ erg 
s~^). Assuming then that the thickness of the cloud AR = Nh/u is less than the distance 
R, we find R < 10^'' cm. A lower limit for the distance of the absorber, assuming it is in the 
form of an outflow, can be obtained equating the velocity of the absorbing material to the 
escape velocity at a given radius R from the central black hole; the derived distance is then 
R ^ lOORg. A constraint on the size of this cloud can be placed assuming 20 ksec as the 
characteristic variability timescale, when our line of sight intercepts the absorbing cloud; this 
gives us AR ^ 6 x 10 cm, w hich corresponds to ~ 10 Rg for a black hole mass of 5 X IO^Mq 
( jWandel fc Mushotzkylll986l ). We can then infer a density of ~ 10^ cm~^ and assuming a 
spherical geometry for the cloud, a mass of about 10^^ g. These one order of magnitude 
estimates for the mass and velocity correspond to a kinetic energy i^km ~ 5 x 10^^ erg and 
using the 20 ksec as a characteristic timescale to a power of ~ 2.5 x 10^^ erg s~^. This value 
corresponds to ~10% of the 2-10 keV X-ray lumino sity and is thus is agree ment with a ra- 



diatio n driven wind model (jProga fc KallmanI 120041 ) or with an aborted jet (jGhisellini et al. 
2004( 1. 



In conclusion, this deep XMM-Newton observation revealed that the soft X-ray emission 
of MCG-5-23-16 can be ascribed to material photoionized by the AGN, likely to be located 
outside the sub-pc scale of the absorber and perhaps coincident with the NLR. We confirm the 
presence of a iron Ka emission line complex composed by a narrow and a broad relativistic 
component. The inclination derived from the diskline profi le ji ~ 40°) is in agreement 
with the orientationally based Unification Scheme of AGN (lAntonucci I Il993l ). the X-ray 
classification of MCG-5-23-16 as Compton thin Seyfert 2 (i.e. intermediate between a type 
1 AGN and a Compton thick type 2) and the optical classification as a Seyfert 1.9. Finally 
we detected a sporadic Fe K absorption feature which could be a signature of a variable high 
velocity outflow. This detection adds one more example to the increasing sample of AGN 
where relativistic outflows have been revealed in the X-ray band. The growing evidence 
of high velocity outflows in AGN indicates that they may play an important role in the 
energetics of AGN central engines. 



We would like to thank the anonymous referee for his/her useful comments that have 
improved this paper. This paper has made use of observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions directly funded by ESA Member 
States and the USA (NASA). Support for this work was provided by the National Aero- 
nautics and Space Administration through Chandra Award Number G05-6146Z issued by 
the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysi- 



-26- 



cal Observatory for and on behalf of the National Aeronautics Space Administration under 
contract NAS8-03060. 

REFERENCES 

Antonucci, R. 1993, ARA&A, 31, 473 

Balestra, I., Bianchi, S., & Matt, G. 2004, A&A, 415, 437 

Bautista, M. A., & Kallman, T. R. 2001, ApJS, 134, 139 

Bianchi, S., Matt, G., Balestra, I., Guainazzi, M., k Perola, G. C. 2004, A&A, 422, 65 
Bianchi, S., Miniutti, G., Fabian, A. C., & Iwasawa, K. 2005, MNRAS, 360, 380 
Bianchi, S., Guainazzi, M., & Chiaberge, M. 2006, A&A, 448, 499 

Brinkman, A. C., Kaastra, J. S., van der Meer, R. L. J., Kinkhabwala, A., Behar, E., Kahn, 
S. M., Paerels, F. B. S., & Sako, M. 2002, A&A, 396, 761 

Cappi, M. 2006, Astronomische Nachrichten, 327, 1012 

Chartas, G., Brandt, W. N., Gallagher, S. C., & Garmire, G. R 2002, ApJ, 579, 169 
Chartas, G., Brandt, W. N., & Gallagher, S. C. 2003, ApJ, 595, 85 
Cash, W. 1979, ApJ, 228, 939 

Dadina, M., Cappi, M., Malaguti, G., Ponti, G., & de Rosa, A. 2005, A&A, 442, 461 
den Herder, J. W., et al. 2001, A&A, 365, L7 

Dewangan, G. C, Griffiths, R. E., & Schurch, N. J. 2003, ApJ, 592, 52 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Fabian, A. C, Rees, M. J., Stella, L., & White, N. E. 1989, MNRAS, 238, 729 

Fabian, A. C, Iwasawa, K., Reynolds, C. S., & Young, A. J. 2000, PASP, 112, 1145 

Ferruit, P., Wilson, A. S., & Mulchaey J. 2000, ApJS, 128, 139 

Ghisellini, G., Haardt, F., & Matt, G. 2004, A&A, 413, 535 

Guainazzi, M., Bianchi, S., & Dovciak, M. 2006, Astronomische Nachrichten, 327, 1032 



Guainazzi, M., & Bianchi, S. 2007, MNRAS, 374, 1290 

Kallman, T. R., Liedahl, D., Osterheld, A., Goldstein, W., & Kahn, S. 1996, ApJ, 465, 994 
King, A. R., & Pounds, K. A. 2003, MNRAS, 345, 657 
Kinkhabwala, A., et al. 2002, ApJ, 575, 732 

Iwasawa, K., Wilson, A. S., Fabian, A. C, & Young, A. J. 2003, MNRAS, 345, 369 
Liedahl, D. A. 1999, LNP Vol. 520: X-Ray Spectroscopy in Astrophysics, 520, 189 
Liedahl, D. A., & Paerels, F. 1996, ApJ, 468, L33 

Longinotti, A. L., Sim, S. A., Nandra, K., & Cappi, M. 2007, MNRAS, 374, 237 
Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, 273, 837 

Markert, T. H., Canizares, C. R., Dewey, D., McGuirk, M., Pak, C. S., & Schattenburg, 
M. L. 1994, Proc. SPIE, 2280, 168 

Markowitz, A., Reeves, J. N., & Braito, V. 2006, ApJ, 646, 783 

Matt, G. 2002, MNRAS, 337, 147 

Matt, G., Bianchi, S., D'Ammando, F., & Martocchia, A. 2004, A&A, 421, 473 

Mattson, B. J., & Weaver, K. A. 2004, ApJ, 601, 771 

McKernan, B., Yaqoob, T., & Reynolds, C. S. 2005, MNRAS, 361, 1337 

McKernan, B., Yaqoob, T., & Reynolds, C. S. 2004, ApJ, 617, 232 

Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197 

Miniutti, G., ct al. 2007, PASJ, 59, 315 

Nandra, K., George, I. M., Mushotzky R. F., Turner, T. J., & Yaqoob, T. 1997, ApJ, 477, 
602 

Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1999, ApJ, 523, 
L17 

Nandra, K., O'Neill, P. M., George, 1. M., Reeves, J. N., & Turner, T. J. 2006, Astronomische 
Nachrich 



-28- 

Perola, G. C, Matt, G., Cappi, M., Fiore, F., Guainazzi, M., Maraschi, L., Petrucci, P. O., 
& Piro, L. 2002, A&A, 389, 802 

Ogle, P. M., Brookings, T., Canizares, C. R., Lee, J. C, & Marshall, H. L. 2003, A&A, 402, 
849 

Porquet, D., & Dubau, J. 2000, A&AS, 143, 495 

Porquet, D., Reeves, J. N., Uttley, P., & Turner, T. J. 2004, A&A, 427, 101 

Protassov, R., van Dyk, D. A., Connors, A., Kashyap, V. L., & Siemiginowska, A. 2002, 
ApJ, 571, 545 

Papadakis, I. E., Brinkmann, W., Page, M. J., Mc Hardy, I., & Uttley, P. 2007, A&A, 461, 
931 

Pounds, K. A., Reeves, J. N., King, A. R., & Page, K. L. 2004, MNRAS, 350, 10 

Pounds, K. A., Reeves, J. N., King, A. R., Page, K. L., O'Brien, P. T., & Turner, M. J. L. 
2003, MNRAS, 345, 705 

Pounds, K. A., & Page, K. L. 2005, MNRAS, 360, 1123 

Proga, D., & Kallman, T. R. 2004, ApJ, 616, 688 

Proga, D., Stone, J. M., & Kallman, T. R. 2000, ApJ, 543, 686 

Reynolds, C. S., & Nowak, M. A. 2003, Phys. Rep., 377, 389 

Reeves, J. N., O'Brien, P. T., & Ward, M. J. 2003, ApJ, 593, L65 

Reeves, J. N., Nandra, K., George, I. M., Pounds, K. A., Turner, T. J., & Yaqoob, T. 2004, 
ApJ, 602, 648 

Reeves, J. N., et al. 2007, PASJ, 59, 301 

Risaliti, G. 2002, A&A, 386, 379 

Risaliti, G., Elvis, M., & Nicastro, F. 2002, ApJ, 571, 234 

Risaliti, G., Bianchi, S., Matt, G., Baldi, A., Elvis, M., Fabbiano, G., & Zezas, A. 2005, 
ApJ, 630, L129 

Sako, M., Kahn, S. M., Paerels, F., & Liedahl, D. A. 2000, ApJ, 543, L115 



-29- 



Sambruna, R. M., Netzer, H., Kaspi, S., Brandt, W. N., Chartas, G., Garmire, G. P., Nousek, 
J. A., & Weaver, K. A. 2001, ApJ, 546, L13 

Schurch, N. J., Warwick, R. S., Griffiths, R. E., & Sembay, S. 2003, MNRAS, 345, 423 

Tanaka, Y., et al. 1995, Nature, 375, 659 

Turner, T. J., Kraemer, S. B., George, I. M., Reeves, J. N., & Bottorff, M. C. 2005, ApJ, 
618, 155 

Vaughan, S., & Fabian, A. C. 2004, MNRAS, 348, 1415 
Yaqoob, T., & Padmanabhan, U. 2004, ApJ, 604, 63 
Yaqoob, T., & Serlemitsos, P. 2005, ApJ, 623, 112 

Young, A. J., Lee, J. C., Fabian, A. C., Reynolds, C. S., Gibson, R. R., & Canizares, C. R. 
2005, ApJ, 631, 733 

Weaver, K. A., Krolik, J. H., & Pier, E. A. 1998, ApJ, 498, 213 

Weaver, K. A., Yaqoob, T., Mushotzky, R. F., Nousek, J., Hayashi, I., & Koyama, K. 1997, 
ApJ, 474, 675 

Weaver, K. A. 2001, in ASP Conf. Ser. 249, The Central Kiloparsec of Starbursts and AGN: 
The La Palma Connection, ed. J. H. Knapen, J. E. Beckman, I. Shlosman, & T. J. 
Mahoney (San Francisco: ASP), 389 

Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 

Wandel, A., & Mushotzky R. F. 1986, ApJ, 306, L61 



This preprint was prepared with the A AS IM^jX macros v5.2. 



-30- 



T3 
O 

E 



o 

T3 




Rest Energy (keV) 

Fig. 1. — Ratio of the XMM-Newton pn data and the model when fitting an absorbed and 
redshifted power law over the 2-10 keV band ignoring the iron emission line energy range 
(5.5 — 7.5 keV). The photon index of the power law has been fixed to 1.8. Two different 
residuals are clearly present: a soft X-ray excess below ~ 2 keV; and the Fe complex at ~ 6.4 
keV. 
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Fig. 2.— The ratio between the XMM-Newton pn (black points), MOSl (red), and M0S2 
(green) data and the the model, plotted between 0.5-1.5 keV. The model consists of the 
baseline absorbed power-law described in the text and an additional soft (F = 3) power law, 
absorbed only by a Galactic column. An emission line is clearly detected at 0.9 keV. 
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Fig. 3.— Left panel: XMM-Newton RGSl and RGS2 spectra in the 0.4-1.3 keV energy 
range. The underlying AGN continuum model is composed of an absorbed power law plus a 
scattered soft power-law component (F = 1.82). Several emission lines are clearly detected. 
The possible identifications of the 5 brightest hnes are also shown. Above ~ 0.9 keV a 
steep rise of the continuum is evident, due to the emergence of the absorbed power- law 
component. Right panel: best fit photoionized plasma model (in the 0.4—1.3 kcV) which 
includes a XSTAR component with an ionization parameter log^ — 1.291q;^3 (^^^ Section 3.2). 
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Fig. 4. — Ratio between the pn data and the absorbed power- law model (F = 1.65) showing 
the iron line profile. The data clearly show a narrow core at 6.4 keV, red and blue wings 
extending from 5.7 keV to 7 keV and a narrow emission feature at ~7.05 keV, which is due 
to Fe K/3. A sharp drop is also present at ~ 7.1 keV due to presence of reflection. 
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Fig. 5. — Data/model residuals for the pn data at the Fe K band. Panel (a) shows the 
data/model ratio when the underlying continuum is composed of an absorbed power law 
plus a reflection component {R = 1.1). Panel (b) shows the residuals left when the narrow 
Fe Ka and Fe K/3 hnes are added to the model. An excess between 5.8 keV and 7 keV is 
still present. Panel (c) shows the residuals when the broad component of the Fe Ka line is 
fitted with a diskline model. The overall fit is now good and only a weak absorption feature 
is left at 8 keV. 
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Fig. 6. — Data/model ratio for the Chandra HETG. The AGN continuum is composed of 
an absorbed power law plus a reflection component (F = 1.82; R = 1.1) A strong narrow 
core at 6.4 keV is observed. When modeled with a single Gaussian, the line is found to have 
EW ~ 80 eV and a width a = 32l|^ eV. 
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Energy (keV) 

Fig. 7. — XMM-Newton spectrum of MCG-5-23-16. The black points denote the EPIC-pn 
data. The red hne is the total best fit model. The AGN continuum model is composed of: an 
absorbed and a scattered power-law components (green lines online version) plus a reflection 
component (light blue line online version). The Fe K line complex (blue line online version) 
is composed of narrow and broad Fe Ka lines and a narrow K(3 line. A soft X-ray emission 
line at ~ 0.9 keV is also shown (magenta line online version). 
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Fig. 8. — Contribution to the the 5.5-8.4 keV for the five pn spectra extracted with a 
time bin of 20 ksec. The Fe Ka line has been parameterized with two Gaussian fines; aU tfie 
parameters of tfie model except tfie primary power-law normalization fiave been kept tied 
togetfier. Tfie only strong deviation in Ax^ is present in tfie tfiird spectrum at 7.7 keV. 
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Fig. 9. — Time resolved spectral analysis. Panel (a): Fe 6.4 keV narrow core intensity (in 
unit of photons cm~^ s~^) versus the time intervals. The dashed lines correspond to the 90% 
confidence level of the normalization of the narrow core measured in the average spectrum. 
Panel (b): same as panel a) for the broad component. Panel (c): Absolute intensity of the 
absorption feature (in unit of photons cm^^ s^^), the line energy has been fixed to the best 
fit value found for the third interval (7.71 keV). Panel (d): MCG-5-23-16 2-10 keV flux (in 
erg cm~^ s~^) . Error bars and upper hmits are at the 90% confidence level. 
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Fig. 10. — Data/model ratio of the pn spectrum of the third 20 ksec interval, with the data 
binned with a constant energy bin of 80 eV. A deep absorption feature is visible at 7.7 keV 
(see text for details) . 
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. 11. — Derived 99%, 90% and 68% confidence contours of the absorption line intensity 
the observer-frame energy. 



-41 - 



1 ' 1 ' 1 ' r 

Model A - log^=3.7, Vg„,=0.1c 




Energy (keV) 
1 ' 1 ' 1 ' r 

Model B - Iog5=3.0, v^„=0 




Energy (keV) 
1 1 1 1 1 ' r 

Model C - log5=1 .5, v„=0.08c 




Energy (keV) 

Fig. 12. — Example of high ionization warm absorber models. Panel(a) shows the best fit 
model: an ionization state characterized by log^ = 3.7~^q'1 erg cm s"^ and a column density 
of ~ 8 X 10^^ cm with an outflow velocity of 0.09 ± 0.01 c. The Fexxvi Lya absorption 
feature shifted to 7.7 keV. fits well the observed dip at this energy. Panel (b): the absorption 
feature at 7.7 keV can be accounted for by assuming a column density of A^h = 10^^ cm"^ 
and an ionization parameter of log^ — 3.0 erg cm s~^. However, the model predicts stronger 
absorption due to the Fe xxv ls-2p absorption line at 6.7 kcV. which is not detected. The 
model in panel (c) shows that lower ionization K/3 absorption line without strong Ka from the 
same species, can be obtained assuming a column density of A'^h = 10^^ cm~^, an ionization 
parameter of log^ = 1.5 erg cm s~^. In this scenario a blueshift of ~ 0.08c would still be 
required to model the absorption line at 7.7keV in the spectrum and the column density 
required to model the EW of the K/3 absorption feature is A^h >> lO^^cm"^; which would 
introduce too much continuum absorption below 6 keV, inconsistent with the pn data. 



-42 - 



Table 1. Log of the observations and exposure times. 



Mission Instrument T(total) (ks) T(„et) (ks) Tstart Tstop 



XMM-Newton 


PN 


131.5 


96.2 


08/12/2005 21:11:28 


10/12/2005 09:17:08 


XMM-Newton 


MOSl 


131.5 


101.6 


08/12/2005 20:42:51 


10/12/2005 09:16:48 


XMM-Newton 


MOS2 


131.5 


102.8 


08/12/2005 20:43:51 


10/12/2005 09:16:53 


XMM-Newton 


RGSl 


131.6 


97.2 


08/12/2005 20:41:37 


10/12/2005 09:18:03 


XMM-Newton 


RGS2 


131.6 


97.2 


08/12/2005 20:41:42 


10/12/2005 09:18:03 


Chandra 


ACIS-S HETG 


30 




08/12/2005 17:41:30 


09/12/2005 02:33:59 


Chandra 


ACIS-S HETG 


20 




09/12/2005 20:52:11 


10/12/2005 03:00:10 



Note. — For XMM-Newton the exposure values reported are total and net exposure time after filtering for high-background 
time intervals. 
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Table 2. XMM-Newton RGSl and RGS2: best fit emission lines required fitting with a 

photoionization model. 



Energy 


Flux 


ID 


EW 


AC 


ELab 


(keV) 
(1) 


(10~^ ph cm"^ s"^) 
(2) 


(3) 


(eV) 
(4) 


(5) 


(keV) 
(6) 




0.68l°:i 


Nvil Lja 


10 ±5 


14.3 


0.500 


U.OU^_Q QQ2 


2.30^}:^^ 


VII Hea 




22.5 


0.561 (f) 
0.569 (i) 
0.574 (r) 






VIII Lya 


14.71?:° 


15.4 


0.654 


n 7Q9+0.007 


0.541°:^° 


Ovil RRC 


1 r. 7+11.0 
^^■'-10.0 


7.0 


> 0.739 


n qnQ+0.020 


0.42l°:i 


Ne IX Hea 


15.9 ±10.0 


10.4 


0.905 (f) 
0.915 (i) 
0.922 (r) 



Note. — The model of underlying AGN continuum has been parametrized 
according to the best fit model obtained with the analysis of the pn spectrum 
(section 3) . The F of the soft power law is tied to the hard power-law component 
and fixed to the best fit value (F = 1.82). The energy of the lines are quoted 
in the rest frame. Fluxes and possible identifications are reported in col. (2) 
and (3) respectively. The EWs are reported in col. (4) and they are calculated 
against the soft power law component. In col. (5) the improvement of fit is 
shown using the C-statistic; the value for the model with no lines is C=455.7 
for 358 PHA bins. In col. (6) we report the theoretical value for the transition 
energies. 
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Table 3. Results of the fit for the mean spectrum and the low and high flux states. 



Parameter Mean High Low 



Continuum T 1.82 ±0.01 1.84 ±0.01 1.81 ± 0.01 

1.49 ±0.01 1.50 ±0.02 1.49 ±0.02 
Flux^ 8.16 8.84 7.39 



Narrow Gaussian E 

EW 

Broad Gaussian E 

a 
EW 

Diskline Rj„ 

i 

EW 



Note. — The line energies are expressed in units of keV, while their- 
widths a and EWs are in eV. The disk radial emissivity has been fixed 
to q — 3. 

^Column density in units of 10^^ cm~^ 
^2-10 keV flux in units of 10~^^ erg cm~^ s~^ 
'^Normalization of the Fe line in unit of 10~^ photons cm~^ s~^ 
^Indicates that the parameter has been kept fixed. 



6.42 ±0.01 6.41 ±0.02 6.40 ± 0.02 



^•^-0.6 



48; 

4.6- 



-12 
-13 
hl.3 
-1.2 



81 



-15 
-15 



6.5^- 



6.4 f 
0.35 ±0.1 

5.9 ±1.5 



6.4 f 
0.351^:11 
611^3^ 



6.4 f 
0.37lg:?l 
75li 
6.311^ 



53: 



hl4 



J-13 

4.6 ±1.2 
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Table 4. Best fit parameters for the absorpion feauture detected at 7.7 keV for the EPIC 

pn and MOS cameras and SUZAKU. 



Mission- Instrument 


Energy 


Flux 


\EW\ 




(keV) 


(lO-^ph cm-2 s-i) 


(eV) 


EPIC-pn'^ 


7.9 ±0.1 


-1.9 ±0.6 




EPIC-pn'' 


7.7 ±0.1 


-3.2 ±0.9 


52± 15 


EPIC-MOS" 


7.4 ± 0.2 


-1.1 ± 1.0 


16± 14 


EPIC-MOS* 


7.4 ±0.2 


-2.2 ± 1.8 


27 ±23 


Suzaku-XISO-1-2-3" 


7.8 ±0.1 


-1.8 ±0.9 


30± 15 


Suzaku-XISO-1-2-3^''^ 


7.65 ±0.4 


-3.5 ±3.0 


52 ±44 



Note. — The parameters of the absorption hne are derived 
adding an inverted gaussian component to the best fit model. 

^Parameters refers to mean spectrum 

'^Parameters are derived for the 20 ksec time interval where the 
line is detected. 

■^The net exposure of this Suzaku spectum is only ~ 6 ksec 



